Non-linear optics with two trapped atoms 
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We show theoretically that two atomic dipoles in a resonator constitute a non-linear medium, 
whose properties can be controlled through the relative position of the atoms inside the cavity and 
the detuning and intensity of the driving laser. We identify the parameter regime where the system 
operates as a parametric amplifier, based on the cascade emission of the collective dipole of the 
atoms, and determine the corresponding spectrum of squeezing of the field at the cavity output. 
This dynamics could be observed as a result of self-organization of laser-cooled atoms in resonators. 



I. INTRODUCTION 



Quantum light sources are an essential element for im- 
plementations of quantum information processing and 
secure telecommunication with quantum optical sys- 
tems P, [H, H El- Experiments have demonstrated sev- 
eral remarkable milestones, thereby opening promising 
perspectives for implementing controlle d g eneration of 
quantum light @, S 0, g I M ED, E EI M ES EH • 
Besides, these studies touch on the fundamental question 
of how macroscopic nonlinear phenomena emerge from 
the dynamics of quantum systems [ItJ ■ A paradigmatic 
example is the optical parametric amplifier 18]. This 
system is usually realized with non-linear crystals in res- 
onators, where the medium response is characterized by 
the dependence of the macroscopic polarization on the 
electric field, and where symmetries of the crystal can 
enhance a certain nonlinear response over others p^ . [2fj| . 
On the other hand, recent theoretical works pointed out 
that a single atom in a suitable setup can constitute an 
efficient non-linear optical medium operating in the quan- 
tum regime [2l|, [22, [H [24[ • A question, which naturally 
emerges from these works, is how these dynamics scale 
up to a macroscopic non-linear medium, and in particu- 
lar what is the microscopic building block exhibiting the 
essential symmetries controlling the order of the medium 
susceptibility. 

In this article we study the non-linear response of a 
medium constituted by two dipoles confined along the 
axis of an optical resonator, and transversally driven by 
a laser, in a configuration like the one depicted in Fig. [1] 
At certain interatomic distances the state of the field 
at the cavity output can exhibit non-classical features. 
In particular, we show that the system response can be 
switched from a parametric amplifier to a Kerr medium, 
just by varying the intensity of the laser field. The valid- 
ity of our analytical predictions are verified by numerical 
simulations which take into account the internal dynam- 
ics of the atoms and their coupling with the quantized 
mode of the resonator. The effect of atomic vibrations 
on the field at the cavity output is estimated using a semi- 
classical model for the atomic motion. Finally, we dis- 
cuss the possibility of obtaining such patterns, operating 
in the quantum regime, as the result of self-organization 



of laser-cooled atoms in the resonator field 129, 3C 




FIG. 1: Two atoms are confined inside a high-finesse optical 
resonator, their dipoles are driven by a laser and couple to a 
mode of the cavity. The quantum state of the field at the cav- 
ity output can be controlled by the interatomic distance inside 
the resonator and the laser intensity and detuning. A detect- 
ing apparatus measures the field at the cavity output. The 
parameters are defined in Sec. [TTJ 



This article is organized as follows. In Sec.|TT]the model 
is introduced and the basic properties are discussed. In 
Sec. lIIII the response of the atomic medium is determined 
as a function of the atomic position inside the resonator, 
when the atoms are driven by a laser. The steady state 
of cavity and atoms is determined for the specific param- 
eter regime in which the system behaves as an optical 
parametric amplifier. In Sec. HVI we consider the effect of 
the center of mass motion on the cavity field by means 
of a semiclassical model. In Sec. [V] we summarize the 
results and discuss some outlooks. The appendices pro- 
vide details of the calculations presented in Sec. |TTT] and 
Sec El 



II. 



THE THEORETICAL MODEL 



We assume two identical atoms of mass M, which are 
confined inside a standing-wave cavity, and localized at 
the position x\ and x-i, respectively, along the cavity axis. 
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We denote by pi and p 2 the corresponding momenta, and 
by H mcc the Hamiltonian determining the dynamics of 
the center of mass in absence of the coupling with the 
electromagnetic field, which has the form 

H mec = A + A + v^x,) (1) 

(2) 

with V(x±,X2) an external potential, which localizes the 
atoms at their equilibrium positions such that they un- 
dergo small vibrations with respect to the cavity-mode 
wavelength. The relevant internal degrees of freedom of 
the atoms are the ground state \g) and the excited state 
\e) of a dipole transition with dipole moment d, which 
is at frequency luq. The dipoles are driven by a trans- 
verse laser field at frequency wj, and couple to a mode 
of the resonator at frequency uj c and wave vector k, as 
displayed in Fig. Q] A detecting apparatus measures the 
field at the cavity output. 



the master equation 

§- t P = -^{H,p}+£ K p + £ 7 p (7) 

EE Cp (8) 

where 

L K p = k (2apa) — a'ap — pa) a) (9) 

is the superoperator which describes noise due to cavity 
decay, and 

C iP = XI \ ( 2rr ;/V T ; - v}<rjP - pa]°j) (10) 

is the superoperator which describes the quantum noise 
due to spontaneous emission. In the superoperator (|10|) 
the term pj accounts for the mechanical effect of the 
spontaneously emitted photon on the atom in Xj , see for 
instance (25j . 



A. Master Equation 

In the reference frame rotating at the laser frequency 
the coherent dynamics of the atoms and cavity mode is 
described by the Hamiltonian H = H mec + H at + H cav + 
H C av-at + H L . The terms 

H at = -HA a)o 3 (3) 

J=l,2 

H cav = -M c a)a (4) 

describe the system dynamics in absence of coupling with 
the electromagnetic field. Here, A = lul — t^o and S c — 
lul — ^>c are the detunings of the laser from the dipole 
and from the cavity frequency, respectively, Oj = \g)j (e| 

the lowering operator of the atom j, o\J its adjoint, and 
a and are the annihilation and creation operators of a 
photon of the cavity mode. The terms 

H L = UQ £ (a] +a j ), (5) 

.7=1,2 

Hat-cav = K ^ g(Xj) (<J<Tj + <j]aj (6) 
.7=1,2 

describe the interaction of the dipoles with the cavity and 
laser fields, respectively, with O the laser Rabi frequency 
and g(xj) the cavity vacuum coupling strength at Xj, 
with g(xj) = gcos(kxj). In Eq. §5§ the laser wave vector 
is orthogonal to the cavity axis. 

Coupling to the external environment gives rise to dis- 
sipation and decoherence, which is described by sponta- 
neous emission of the excited state at rate 7 and by cavity 
decay at rate re. The dynamics of the density matrix p 
of the cavity and atomic degrees of freedom is given by 



B. Multi-photon processes and atomic patterns 

It is instructive to consider the dynamics in terms of 
the collective states of the dipole. We denote by |+) 
and |— ) the Dicke symmetric and antisymmetric states, 
respectively, with |±) = (\eg) ±\ge)) / V2, and rewrite the 
interaction of the atoms with laser and cavity mode in 
terms of the operators 

S ± = (a 1 ±a 2 )/V2. (11) 

In this representation, the laser-atom interaction, Eq. |5|). 
is rewritten as 

H L = hV2QS+ + H.c. (12) 

while the atom-cavity interaction term, Eq. ([6]), can be 
decomposed as H at - ca v — H + + with 

H± = hg±(xi, X2) (aS± + a)S±) (13) 

and 

g±{xi,x 2 ) = —/= (cos(fcxi) ± cos(fca;2)) . (14) 
v2 

This decomposition highlights the relevant cavity-atom 
dynamics, which depend on the relative atomic posi- 
tion. The term H_ describes the coupling of the cav- 
ity mode with the Dicke anti-symmetric state, and it 
vanishes when the interatomic distance d — x 2 — x\ is 
an integer multiple of the cavity wavelength A = 2n/k. 
We denote the corresponding atomic configuration as " A- 
spaced pattern" . The term H + describes the coupling of 
the cavity mode with the Dicke symmetric state and it 
vanishes when d is an odd multiple of A/2. We denote 
the corresponding atomic configuration as "A/2-spaced 
pattern" Below we discuss the corresponding dynamics 
in detail. 
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1. X-spaced pattern. 

We first consider the case in which the interatomic dis- 
tance is an integer multiple of A. For this configuration, 
at steady state and for large cooperativities, the atoms 
are in the ground state and the cavity mode is in a co- 
herent state whose amplitude is determined by the laser 
intensity [26|, [27} • This behaviour can be understood in 
terms of the coherent buildup of a cavity field, such that 
its phase is opposite to the driving field. As a result, 
the atomic dipole is not excited, even if the cavity mode 
is in a coherent state with a finite number of photons. 
When two or more atoms are present inside the resonator, 
this situation can be achieved when the atoms scatter in 
phase into the cavity modes, i.e., when they are arranged 
in a A-spaced pattern. The coherent scattering processes 
which two atoms undergo are sketched in Fig. [2Ja) in the 
Dicke basis, showing that the antisymmetric state |— ) 
remain always decoupled from the coherent dynamics. 
Here, one identifies the suppression of excitation of the 
atoms at steady state as due to interference between the 
excitation path \gg, n) — > |+, n), driven by the laser, and 
the excitation path \gg,n + 1) — > |+, n), driven by the 
cavity. Figure [2^a) displays also the other higher-order 
processes. In particular, we note the processes which 
lead to the excitation of the state |ee,n) by the absorp- 
tion of two laser photons, followed by emission of pair of 
photons into the cavity. These processes are expected to 
give rise to squeezing of the coherent state of the cav- 
ity field. We note that squeezed-coherent radiation has 
been predicted in the resonance fluorescence of an atomic 
crystal, at wave vectors such that the Bragg condition of 
the atomic crystal is equivalent to the A-spaced pattern 
here discussed (28|. Finally, we note that the formation 
of A-patterns of laser-cooled atoms inside of resonators 
has been predicted as the result of a self-organizing pro- 
cess [1^, [3(| , and features of the field at the cavity out- 
put, associated with their formation, have been measured 
in [3l|, [32j ■ Theoretical works have shown that these pat- 
terns can be also stable in the strong coupling regime, 
under the condition, in which atomic excitation is sup- 
pressed and the cavity field is in a coherent state [III [33[ • 



2. X/2-spaced pattern. 



collective dipole can emit photons in pairs into the cav- 
ity mode. These processes are expected to give rise to 
squeezing of the state of the cavity field. We note that 
squeezed radiation has been predicted in the resonance 
fluorescence of an atomic crystal, at wave vectors such 
that the Bragg condition of the atomic crystal is equiva- 
lent to the A/2-spaced pattern here discussed [28j , In this 
paper we will investigate the quantum state of the light 
in presence of a high-finesse cavity when the atoms are 
initially in a A/2 spaced pattern, and determine the dy- 
namics resulting from the competition between coherent 
processes and noise, such as cavity decay, spontaneous 
emission, and atomic vibrations at the equilibrium posi- 
tions. 




66,n+ 2) 




We now analyze the case, when the interatomic dis- 
tance is an odd multiple of A/ 2, such that H + = 0. In 
this case the atomic ground state couples via the laser 
to the Dicke symmetric state |+), and via the cavity to 
the antisymmetric state |— ), as depicted in Fig. [2jb). 
Hence, when the laser drives the atoms well below sat- 
uration,the cavity is empty 26]. In fact, in this limit 
the two atoms scatter the laser photons with opposite 
phase into the cavity and the resulting field vanishes due 
to destructive interference. Figure HJb) shows, however, 
that the cavity mode can be pumped by higher-order pro- 
cesses, which excite the state |ee,0). In this regime, the 



FIG. 2: Sketch of the coherent scattering processes between 
the collective states of two atomic dipoles driven by a laser and 
coupled to the cavity mode, when (a) the interatomic distance 
d is an integer multiple of the cavity-mode wavelength A, and 
(b) when d is an odd multiple of A/2. The states |J, n) are 
the Dicke states of the two dipoles |J) at n cavity photons, 
where | J) = \gg), |±) , |ee), and |±> = (\eg) ± \ge))/y/2. The 
arrows labeled by Q. (g±) indicate the transitions driven by 
the laser (the cavity mode). 
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III. NON-LINEAR RESPONSE OF TWO 
TRAPPED ATOMS 

In this section, starting from Eq. (0 we derive the 
equation describing the effective dynamics of the cavity 
mode in the limit of large atom-laser detuning |A| ^> 
g, f2, \6 C \, 7, k. In this analysis we neglect the effect of 
atomic motion, and identify the parameter regime in 
which the system operates as a parametric amplifier. The 
prediction of the analytical model are compared with the 
results of a numerical simulation, which evaluate the cav- 
ity mode state by solving Eq. ([7]). 



A. Effective Hamiltonian 

We derive the effective Hamiltonian H c g for the coher- 
ent cavity dynamics at fourth order in the expansion in 
the small parameters <?/|A|, f2/|A|. In the Hilbert sub- 
space subtended by the states \gg, n), with n the number 
of cavity photons, it has the form 

H cS = (0 - S c ) a^a + (3 (a f + a) 

Xa) a) aa + ^ {a^ 2 + a 2 ^j , (15) 



where 



9+ (xi,x 2 ) +g 2 _ [xi,x 2 ) 



R ^ ( \ 

P = —z-9+{xi,x 2 ) 



(16) 
(17) 



1 2 

X = -^[9 2 +{x ll x 2 ) - g 2 _{x ll x 2 )] (18) 

a = -^■[9+(xi,X2)-g_(x 1 ,x 2 )] (19) 

Here, 9 is the a. c. -Stark shift experienced by the cav- 
ity field due to the interaction with the atoms, the term 
(3 comes from the H + term, Eq. (|13[) . and results from 
the two-photon transitions coupling the photon states 
\n) and \n ± 1), see Fig. [U(a). The amplitude a is the 
strength of the effective nonlinear pumping of the cav- 
ity field which gives rise to a x^ nonlinearity, typical 
of a degenerate parametric amplifier [l8j |. This term 
is the sum of two contributions, which are weighted by 
g + and g_, respectively, and which represent the coher- 
ent sum of the four-photon processes coupling the states 
\gg, n) — > \gg,n± 2) and depicted in Fig. [2ja) and (b). 
Finally, the amplitude x is the a.c-Stark shift associated 
with four-photon processes, where two cavity photon are 
virtually absorbed and then emitted along the transition 
\gg,n) — > \ee,n — 2). This term is present in both pat- 
terns, and gives rise to the x^ nonlinearity typical of a 
Kerr medium. 

The form of Hamiltonian (|15| highlights how the two 
patterns we considered, A- and A/2-spaced, contribute to 
the various nonlinear processes. We first notice that in 



presence of only one atom (when, e.g., g(x 2 ) = 0) the 
terms a and x trivially vanish: these types of nonlineari- 
ties can be clearly generated only when both atoms cou- 
ple to the cavity mode. Then, one observes that the two 
patterns gives rise to different nonlinear dynamics. In the 
A-spaced pattern, for instance, all terms in Eq. (|15|) con- 
tribute to determine the coherent dynamics of the cavity 
mode. While the linear shift 8 can be set to zero by 
properly choosing the detuning <5 C , on the other hand 
the linear term scaling with (3 is dominant, and one rea- 
sonably expects that it will determine the cavity steady 
state. 

When the atoms are distributed in a A/2-spaced pat- 
tern, the linear drive in Hamiltonian (|15p vanishes, i.e., 
(3 = 0, while the only terms which contribute to the co- 
herent dynamics are at fourth order in the perturbative 
expansion. Two possible scenarios can be here identified, 
(i) When the laser drive is much weaker than the cavity 
coupling, fl <C g, then |x| 3> \a\ and the dynamics will be 
basically equivalent to a Kerr medium as in [2l| , whereby 
in our case the Kerr nonlinearity emerges from the in- 
teraction of the cavity field with the collective dipole of 
the atoms, (ii) When the laser drive is much stronger 
than the cavity coupling, Q 3> g, then \\\ <C |a| and the 
dynamics will be essentially equivalent to the one in a 
X^ 2 -* -medium. This is the case on which we focus in the 
rest of this article. 



B. Realization of a x' 2 ' medium 

We now consider Hamiltonian (|15|) when the atoms 
are localized at the antinodes of the cavity modes in a 
A/2-spaced pattern, i.e., when (3 = 0, and when Q g, 
i.e., \a\ ^> \x\. Setting 5 C = (9, the effective coherent 
dynamics of the cavity mode is described by Hamiltonian 
H eS « H', with 



*'=f(°*w) 

and a = a, where now 

4f>V 

a = TT~ 

A 3 



(20) 



(21) 



A master equation for the reduced density matrix g of the 
cavity mode can be derived from Eq. (|7|) , which takes the 
form 

d i 

— g = -- [H , q] + C K g + CyQ 



where superoperator C K is defined in Eq. ([9]), while 
C^g = — (2aga^ — a) ag — ga^a) 



(22) 



describes the damping of the cavity mode via sponta- 
neous emission, with 7' « 7<? 2 /A 2 . 
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When a > k + 7', Eq. (|22|) predicts that the energy 
of the cavity mode increases exponentially as a function 
of time. Clearly, this exponential increase is a good ap- 
proximation only for short times, when the number of 
photons inside the cavity mode still warrants the validity 
of the perturbative expansion, while for longer times the 
dynamics will be determined by competition with other 
processes which we neglected in the derivation. 

When k + 7' > a, a steady state solution exists, and 
the corresponding stationary average photon number is 



n 



a 



st 



(23) 



where k' = 7' + k. In this case, the field quadrature 

X(t) =a(t)e- i?i +a t (t)e^ (24) 

is squeezed for <j> — 7r/4, and its steady-state variance, 
(AX| t ) = {X 2 ) st - (X) 2 st , takes the form 

<AA s 2 t ) = 



a 



(25) 



Hence, in this case the reduction of the noise of the 
quadrature at steady state is such that ^AJTg t ) > |, 
since k! > a. 

We now identify parameter regimes in which these dy- 
namics can be found. Master Equation (f2"2")) has been 
determined by evaluating the coherent processes up to 
fourth order, treating cavity decay at lowest order, and 
spontaneous emission at second order in the perturba- 
tive expansion. In particular, by deriving the superop- 
erators in Eq. ^ and Eq. (|22p we neglected dissipative 
scattering processes at higher order in the expansion in 
Q/\A\,g/\A\. This is valid provided that g 2 /\A\ > k,j 
and when a > 7', which corresponds to the condition 



7 



< 



n 2 



(26) 



where we used Eq. (|21[) . For a dipole transition with 
linewidth 7/27r = 100 kHz, in a cavity with g/2?t — 
2.7MHz, setting Sl/2n = 10 MHz, |A| /2tt = 100 MHz. 
we find \a\/2ir ps 3 kHz and a negligible rate of spon- 
taneous decay. Appreciable squeezing could be observed 
for a cavity decay rate of few kHz, which is a demanding 
experimental condition. We will focus on this parame- 
ter regime and check numerically the correctness of the 
predictions of our analytical model. 



C. Squeezing spectrum at the cavity output 

Assuming that the system is in the regime where k! > 
a, we evaluate the spectrum of squeezing of the field at 
the cavity output, namely [l8j 

/>oo 

S ont (uj) = 2Re / dt e~ ia,t (27) 
Jo 

x «X out (t)Jf out (0)> st - (X ont (t)) st (X out (0)) st ) 



where the subscript St indicates that the averages are 
performed over the steady state density matrix. In 
Eq. (f2"7| X out (t) is the quadrature of the output field, 
defined as 



X out (t) = OoutWe" 1 * + aLt We i0 
with <f> = 7r/4 and where 

flout (*) = \/Ka(t) - a in (t) 



(28) 



(29) 



and di n (t) is the input noise which is delta-correlated, 

(^ai a (t)al n (t')) — S(t — t'). Using the effective model in 

Eq. (|22p we find an analytical expression of the squeezing 
spectrum 



SoutM = 1 - 



Ana 



(30) 



showing that a large reduction of the quadrature fluctu- 
ations below the shot noise limit is achieved at u = 
when k' w a. 

Figure[3]displays the spectrum of squeezing, comparing 
the analytical prediction in Eq. (|30[) with the numerical 
result obtained using Eq. ([7]), hence including the full 
internal dynamics of cavity and atoms, as well as the in- 
coherent processes due to cavity decay and atomic spon- 
taneous emission at all orders, as discussed in App. [A] 
The spectra are evaluated by setting a = k/2, and show 
that for this parameter regime the analytical model pro- 
vides a good description of the dynamics. We note, as 
expected, that spontaneous emission tends to decrease 
the squeezing at the cavity output. Figure [4] displays the 
spectra of squeezing for a larger value of the cavity cou- 
pling strentgh. Discrepancies between the analytical and 
the numerical model arise from the contribution of the 
Kerr non-linearity in Eq. (|15p . which is not negligible for 
this parameter regime, since the laser Rabi frequency Q 
and the cavity coupling strength g are of the same order 
of magnitude. 

Figure [5] displays the value of the squeezing spectrum 
at to = as a function of the cavity decay rate k. The 
spectrum is plotted for k > a, when the analytical model 
described by Eqs. (|22|) allows for a steady state solution, 
and it clearly shows that squeezing at the cavity output 
is very sensitive to variations of k. On the other hand, 
the dependence on the atomic linewidth 7 is compara- 
tively weak, as one can see from Fig. [6] The discrepancy 
between numerical and analytical model at lower values 
of 7 is due to the contribution of incoherent scattering 
processes at higher order, which are accounted for in the 
numerics and give rise to a very narrow peak at to = in 
S(lj). This feature however does not appear for shorter 
integration times, corresponding to the limit of validity 
of our perturbative treatment. 

Figures [7^ a) and (b) display the spectrum of squeezing 
at lu — and the corresponding variance of the maxi- 
mally squeezed quadrature of the cavity field as a func- 
tions of k' = k + 7'. In Fig. [T^a) the upper curves are 
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FIG. 3: Squeezing spectrum of the field at the cavity out- 
put, when the atoms are in a A/2-pattern. The dashed 
lines correspond to the spectrum evaluated analytically from 
Eq. (|30|l . the solid lines to the spectrum found from the nu- 
merical evaluation of the steady state of Eq. 0, see App. lAl 
The frequency is in units of reo = K- The parameters are 
A = -1.25 x 10 5 re , Q = 1.25 x 10 4 re , g = 1.25 x 10 3 re 
and S c = — 24ko (<5 c is chosen so to compensate all a.c.-Stark 
shifts). For the choice of these parameters, a — k/2. The 
lower and upper curves have been evaluated for 7 = and 
7' = reo/2 (7 = 10 4 kq), respectively. 




FIG. 5: Value of the squeezing spectrum S ut(u>) at u — 0, 
S'out(O), as a function of the cavity decay rate re in units 
of kq. The dashed lines correspond to the value predicted 
from Eq. (|30|l . the solid lines to the numerical result found 
from Eq. ©. The parameters are g = 1.25 x 10 3 «o, A = 
-1.25 x 10 5 re , n = 1.25 x 10 4 re , and 5 C = -24fto. For these 
parameters a = kq/2. The lower and upper curves have been 
evaluated for 7 = and 7' = reo/2 (7 = 10 4 reo), respectively. 
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FIG. 6: S'out(O) as a function of the atomic spontaneous emis- 
sion rate 7 in units of reo. The spectra are plotted for two 
values of the cavity decay rate re = reo and re = /to/2. The 
other parameters are as in Fig. 



FIG. 4: Same as in Fig. [3] where now re = lOOreo, g = 
1.25 x 10 4 reo, and 8 C = —24 x 10 2 re ( ). For the choice of 
these parameters, a = re/2. The discrepancy between ana- 
lytical and numerical results is due to the contribution of the 
Kerr-nonlinearity, which is not accounted for in the analytical 
model. 



obtained for k = 7' = k'/2, the lower curves correspond 
to 7' = 0, k' = K. Figure \T[b) shows that the vari- 
ance of the quadrature is the same both for 7=0 and 
7' = k, showing that spontaneous emission in this regime 
only dissipates the squeezed field along other channels, as 
predicted from the analytical model of Eq. 



IV. EFFECT OF THE ATOMIC MOTION 

So far we have studied the dynamics of the cavity mode 
neglecting the atomic kinetic energy on the cavity-mode 
dynamics. In this section we study the effect of fluctua- 
tions in the atomic positions, when the system operates 
as an optical parameteric amplifier. We assume that the 
atoms are confined by an external potential, which lo- 
calize them at the antinodes of the cavity standing wave 
in a A/2-spaced pattern, in the regime in which the me- 
chanical effects of the cavity field on the atomic motion 
can be neglected. This situation could be realized ex- 
perimentally with the technology developed for instance 

in MM, M M. 
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FIG. 7: (a) S'out(O) as a function of the total effective dis- 
sipation rate k' — k + 7' in unit of «o and (b) correspond- 
ing variance of the squeezed quadrature of the cavity field. 
The numerical results are displayed for k — k', 7 = (cir- 
cles) and k = 7' = k'/2 (crosses). The dashed lines are 
obtained from the analytical model. The other parameters 
are g = 1.25 x 10 3 k , S c = -24k , A = -1.25 x 10 5 k , 
D. = 1.25 x 10 4 fto, and a = Ko/2. 



Denoting by Xj the atomic equilibrium positions, and 



by Qj 



• Xj the displacements, we write the external 



potential for small vibrations as 



1 



V{X U X 2 ) = -MS (qi + q A 2 ) 



(31) 



where v is the trapping frequency. The Hciscnbcrg- 
Langevin equation of motion for the atomic displacement 
q.j is given by [38| 



(32) 



where £(f) is the quantum Langevin force, associated 
with the spontaneous emission and the cavity decay pro- 
cesses, and 



OX j 



(33) 



is the mechanical force operator arising from the spatial 
gradient of the atom-cavity interaction over the atomic 



wave packet. These equations have to be solved to- 
gether with the Heisenberg-Langevin equations for the 
field, which depend on the atomic motion through the 
functions cos kxj . We assume that the atoms are well lo- 
calized at the antinodes of the cavity mode, namely that 

5q = \J (<lj) *C A, and make a perturbative expansion in 
the small parameter kSq. At second order, the equations 
for the fields read 



d(f) = -iaa f (t)- [k +i(9 - 8 c )]a(t) + r)(t) 



+iy [(«? + «§) (aa\t) + ea{t)) 
a't(t) = iaa(f)- [k' - i (6> - <5 C )] a f (f) + 77 (f) 



(34) 

<&) a] 

(35) 



-iy M + ii) («»(*) + ^ f (*)) + (9? - ii) 0} 

with a defined in Eqs. ((21]). = gQ/A, 8 = 2g 2 A, and 
r](t) is the quantum Langevin term, 77(f) — ^/2~K,a m (t) + 
%/27'afj* (f ) . Here af* (f ) is the input noise term associated 
with atomic spontaneous emission, which satisfy the re- 

lation (<(t)o£ t (t'))=*(*-* / ). 

Even when the atoms are well localized around the 
antinodes of the cavity mode, the systematic solution of 
these coupled equations is rather complex. Here, we as- 
sume that the external potential provides a steep con- 
finement, such that the effect of the coupling with the 
cavity mode can be neglected in Eq. (f3"2"|) . In this limit 
the solution of Eq. (f32|) reads 



OS it) 



qf 1 ' cos(^f + 



(36) 



where q^ 



and 



, j)j are determined by the initial condi- 

tions. When the trap frequency is much larger than the 
effective rates which determine the evolution of the field, 
v 3> a, k', we can derive a secular equation for the cavity 
field by averaging the equations for the cavity variables 
over a period T — 2n/v [3§|. We insert Eq. (|36jl into 
the equations for the field variables, Eq. (j3~4"|) - ([33]) . and 
integrate them over the period T. With this procedure 
we find equations for the operators a(f), 77(f), defined as 



a(t) 



t+T 



dra(r), 77(f) = 



t+T 



dTT)(T) (37) 



Here, the new noise operators satisfy the equation 
(fj(t)tf (t')) ~ 2n'8{t — f'), where the S-\ike correlation is 
to be interpreted for the coarse-grained time scale. The 
corresponding Heisenberg-Langevin equations read 

a(t) = -iaa f (f)- K + i(e-S c J a(f) + 77(f) (38) 
a(f)t = iaa(t)- Ik - i (o - S c )] tf(t)+tf(t) (39) 



(40) 



while their derivation is discussed in App. [B] Here, 

a = a{l-k 2 q 2 /2) 
6 = 6{l~k 2 q 2 /2). 
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and we have assumed that the oscillation amplitudes of 
the two atoms are equal, gj ^ = q^ = q. The motion- 
induced a. c. -Stark shift can be compensated by prop- 
erly tuning the laser frequency, 5 C = 0(1 — k 2 q 2 /2), and 
Eqs. (|38 ]) -(|39] 1 become 



~a(t) 
5(*) 



= -\aa\t) - n'a(t) + fj(t) 
= iaa(t) - Ka\t) +fj ] (t). 



(41) 



Correspondingly, at lowest order in kqj the spectrum of 
squeezing is 



SoutM = 1 - 



Ana 



(k' + a) + lu 2 



2/2 9 
(X — K — LU 



k 2 q 2 



(K' + a) 2 + LJ 2 



(42) 



where the term proportional to k 2 q 2 is the correction 



to Eq. (|30p due to small vibrations of the atoms at the 
equilibrium positions. Small fluctuations hence reduce 
the bandwidth of frequencies where the light is squeezed. 
The corresponding spectrum, Eq. (|42p , is displayed in 
Fig. [8] for kq = 0.3 and compared to the one of Eq. (|30"|) . 
where atomic motion is neglected, showing that the mod- 
ification of the spectrum of squeezing due to the motion 
is very small. 




the intensity of the laser the system can operate as Kerr 
medium or as optical parametric amplifier, whereby the 
nonlinear response emerges from the collective excita- 
tions of the atomic dipoles. We have studied in detail the 
case in which the system operates as an optical paramet- 
ric amplifier, and investigated the squeezing of the field at 
the cavity output considering the effects of atomic vibra- 
tions, when the atoms are confined inside the resonator 
at the equilibrium positions of a steep external potential, 
in a situation which can be experimentally realized for 
instance in H [M HE [53 . 



A natural question, emerging from recent studies on 
selforganization of laser-cooled atoms in resonators [29|, 
l30l . [3ll [3^ , is whether in absence of an external potential 
trapping the atoms, the A/2-spaced pattern can be sus- 
tained by the mechanical forces of the potential generated 
by the scattered field. In [33[ a semiclassical and numeri- 
cal analysis showed that this configuration is expected to 
be stable for choices of the parameters, which are consis- 
tent with the operational regime in which squeezed light 
can be observed. In this case, one would hence have a 
selforganized pattern, which sustains and is sustained by 
non-classical light. 

The results of this work provide an example of how 
non-linearities emerge from the microscopic dynamics of 
few simple quantum systems. In this respect, two atoms 
in a resonator can be considered the most basic realiza- 
tion of a non-linear crystal, with however limited effi- 
ciencies. We conjecture that by scaling up the number of 
atoms collective effects can enhance the nonlinear prop- 
erties, thus improving the system response. Another in- 
teresting question is how the system dynamics are mod- 
ified when t he q uantum nature of the atomic motion is 
relevant [iol |4l| . and in particular how the correlation 
functions of the output field are affected by the quantum 
properties of the medium. This study requires an anal- 
ysis of the spectrum of resonance fluorescence as in (42j, 
which systematically accounts for the quantum state of 
the atomic motion, and it will be object of future inves- 
tigations. 



FIG. 8: Spectrum of squeezing of the field at the cavity output 
as a function of u in units of Ko, for the same parameters as in 
Fig.0and 7 = 0. The solid curve corresponds to the spectrum 
of Eq. (|4ip for kq = 0.3. The dashed line corresponds to the 
spectrum of Eq. (pO) when atomic vibrations are neglected. 
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APPENDIX A: EVALUATION OF THE 
SQUEEZING SPECTRUM 

Using Eq (|29p . we rewrite the squeezing spectrum in 
Eq. (22) as 



V 



i6/2 ia/2 
-ia/2 -i0/2 



and 



B 



-i/3 



S'out(w) = 1 + 4«;Re / cite 
Jo 



(a(*),a(0))ste 



3 -2i7r/4 



We indicate with 



+( fl f(0) I at(t)) st e 2i ^ + ( fl t(t), ft ( )) St + 



;ot(o),o(t))st 

(Al) 



/(*) = 5; 



t+T 



dr/(r) 



(B5) 



(B6) 



(B7) 



where 



( a t( ), a t(t)) st = ( a t( ) a t(t)) st - (at(0)) S t(a t (t))st. 

Equation (|A1[) can be expressed in terms of averages 
performed over a density matrix by means of the rela- 
tion (A(t)A(0)) St = Tr {Ae a Ap st } and (A(0)A(t)) st = 
Tr {Ae Ct pstA} , where A is a generic operator, L is the 
Liouvillian defined in Eq. (SJ) setting H mec = 0, and pgt 
is the steady state density matrix satisfying the relation 
£pst = 0. Therefore the spectrum of squeezing can be 
rewritten as 



sS( W ) = l-4«Re 



7r5(o;)Tr{xWp st y 



+ Tr {x^ (C - ilo)- 1 (a Pst e- ie + p st ah w ) }JA2) 



The numerical results in Sec. IIII Cl are based on the eval- 
uation of the spectrum of squeezing, as calculated from 
Eq. (|Al2|) using the Liouvillian of Eq. ([7}. 



APPENDIX B: DERIVATION OF THE SECULAR 
EQUATIONS FOR FAST VIBRATING ATOMS 

After inserting Eq. ||36J) into the Eqs. Eqs. (|M ]) -([55 |) . we 
obtain 

A(t) = MA(t) + N(t) 



3 = 1,2 



where 



a(t) 
ot(t) 



(Bl) 



(B2) 



N(t) 



vHt)J> 



(B4) 



the time average of a variable f(t) over a period of oscil- 
lation T = 2~n jv of the atomic motion. Since 

|/» = i[/(* + T)-/( t )] = |, 

we find 



Mi(t) + N(t) + k 2 q 2 VA{t) 

rt+T 



+ k2 ^ J2 Y I dTCOs(2vT + 2<l> j )A{T) 



3=1,2 



(B8) 



where we have used the relation cos 2 (y) = h [1 + cos(2y)] 
and we have assumed that the two atoms have the same 



energy, such that q\ 



q . We now identify the 



conditions under which we can neglect the second line of 
Eq. (|B8j) . Integrating by part the second line of Eq. (|B8|) 
an using Eqs. (|B1[) and (|B8[) we obtain 



MA(t) + N(t) + k 2 q 2 VA(t) 
+k 2 q 2 C{t) + k A ^D(t) + k A <fE (B9) 



where 
C(t) - 



^^|siri(2i4 + 2c/)j) MA(t)+N(t) 



i r t+T 1 

-- y dr sin(2j/T + 2&) [Mi(r) + N(t)] \ 



D(t) 



8vT 



E 

jj' 



t+T 



ch 



E = 



x [sin(2i4 + 2(f) j) cos(2^r + 2<j)ji) 
-2sin(2^r + 2<f>j) cos 2 (t/r + <pj>)] A(t), 
VB 

— sin 2(02 • 
8^ 



(BIO) 



The terms k 2 q 2 C(t)+k 4 q 4 D(t) are negligible with respect 
to fc 2 g 2 y when |0| k//8i/ < |a| /2 and |6»(6» - S c )\ /8u < 
|a| /2 which reduce to 



v > 
v > 



il 
|A|' 

fc 2 g 2 



g 2 A 



n 2 



(BH) 
(B12) 
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when \a\ and k' are of the same order of magnitude 
and S c = 0(1 - k 2 q 2 /2), see Eq. (gOJ) - The term k 4 q 4 E 
in Eq. (|B9j) can be neglected when k 4 q 4 \6f3\ /IQv <C 
k 2 q 2 \a\ /2, that is 



v > 



k 2 q 2 



(B13) 



When conditions (|B11|) - (|B13|) are satisfied we approxi- 
mate Eq. (|B8|) with 



9 r 



Of 



A(t) = M A(t) + N(i) + k q VA(t). 



(B14) 



which then leads to Eqs. (|38[) and Finally we show 

that the averaged noise operators, fj(t) and ?7^(i), which 
appear in the term N(t), are delta correlated. The only 



non- vanishing correlation function is 

{WW = *"/ dr'6(T-r') 

(t' + T-t) for t' <t <t' + T 

2k' 



„ 2 (i + T - t') for t' - T < t < t' 
for t > t' + T or t < t' - T 



(B15) 



which is not zero only if the two integration intervals 
[t, t + T] and [t 1 , t' + T] have finite overlap. Therefore if 
f(t) varies slowly over the time T, so that f(t) ~ fit), 
then one has 



(B16) 



dtf(t)(fj(t)tf(t'))^2K'f(t'), 

) 

that is (fj(t)tf(t')) ~ 2k'<J(* - t'). 
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